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Abstract 

We present a shell model calculation for the beta decay of ^^C to the ^^N ground- 
state, treating the relevant nuclear states as two Op-holes in an ^^O core. Employing 
the universal low-momentum nucleon-nucleon potential V|ow-k only, one finds that the 
Gamow- Teller matrix element is too large to describe the known (very long) lifetime of 
^^C. As a novel approach to the problem, we invoke the chiral three-nucleon force (3NF) at 
leading order and derive from it a density-dependent in-medium NN interaction. Including 
this effective in-medium NN interaction, the Gamow- Teller matrix element vanishes for 
a nuclear density close to that of saturated nuclear matter, po = 0.16fm~^. The genuine 
short-range part of the three-nucleon interaction plays a particularly important role in this 
context, since the medium modifications to the pion propagator and pion-nucleon vertex 
(due to the long-range 3NF) tend to cancel out in the relevant observable. We discuss 
also uncertainties related to the off-shell extrapolation of the in-medium NN interaction. 
Using the off-shell behavior of Viow-k as a guide, we find that these uncertainties are 
rather small. 



PACS: 21.30.Fe, 21.60.Cs, 23.40.-s 
Keywords: Effective field theory at finite density, chiral three-nucleon force, shell-model calculation of 
^^C beta-decay. 

1 Introduction 

The anomalously long beta decay lifetime of ^^C, which makes possible the radiocarbon dating 
method, has long been a challenge to nuclear structure theory. The transition from the Jf = 0^, 
Tj = 1 ground state of ^^C to the JJ = l"*", Tf = ground state of ^^N is of the allowed Gamow- 
Teller type, yet the known lifetime of ~ 5730 ±30 years is nearly six orders of magnitude longer 
than would be expected from typical allowed transitions in p-shell nuclei [1, 2]. The associated 
Gamow- Teller (GT) transition matrix element for the ^^C decay must therefore be accidently 
small, on the order of 2 x 10~^, which makes this transition a sensitive test for both nuclear 
interaction models and nuclear many-body methods. 

The earliest studies to give insight into this problem were based on phenomenological models 
of the residual nuclear interaction. If only central and spin-orbit forces are included, it can be 
shown [3] that it is impossible to achieve a vanishing GT matrix element when the model space 
is restricted to Op^^ configurations consisting of two Op holes in an ^^O core. Later, Jancovici 
and Talmi discovered [4] that this problem can be overcome with the addition of a tensor force 
component in the residual interaction. When realistic nuclear forces based on meson exchange 
were developed in the 1960's and applied to nuclear many-body problems through the G-matrix 
formalism, Zamick found [5] the decay rate to be very sensitive to the Opi/2 — 0p3/2 splitting, 

^Work supported in part by BMBF, GSI and by the DFG cluster of excellence: Origin and Structure of the 
Universe. 



1 



with overall unsatisfactory results when the experimental value of 6.3 MeV is used. A more 
recent calculation [6] of the ^''C beta decay has been performed in a large-basis no-core shell 
model study performed with the Argonne V8' interaction [7] . Although it was shown that the 
inclusion of additional configurations up to Ghcu lead to a suppression of the GT matrix element, 
the results were not yet converged at this order. 

Very recently, it has been suggested [8] that the ^"^C beta decay transition matrix element 
should be particularly sensitive to the density-dependence of the nuclear interaction. The 
study in ref. [8] used a medium-dependent one-boson-exchange nuclear interaction modeled 
with Brown- Rho scahng [9, 10], in which the masses of the bosons (except the pion) decrease in 
a nuclear medium due to either normal hadronic many-body effects or the partial restoration 
of chiral symmetry at finite density. In the present work, we examine the role of density- 
dependent corrections to the nuclear interaction due to the leading-order chiral three-nucleon 
force (3NF) at one-loop order. For the density-independent two-body part of the interaction 
we use the low-momentum interaction Viow-k derived from the Idaho N3L0 chiral potential 
[12, 13] for cutoffs between Aiow-k = 2.1 fm~^ and 2.3 fm~^. For densities in the neighborhood 
of saturated nuclear matter (po = 0.16 fm~^, we find a large suppression of the GT matrix 
element that is almost entirely due to the short-range component of the chiral 3NF. There are 
three contributions arising from the long-range two-pion exchange component of the chiral 3NF, 
but the largest two terms approximately cancel. Moreover, we find that the density-dependent 
corrections resulting from the medium-range 3NF are naturally small. We also calculate the 
Gamow- Teller strengths from the ground state of ^^N to the excited states of ^^C and find that 
they are in satisfactory agreement with recent experimental data [11], although generally they 
are much less sensitive to the density dependence of the nuclear interaction than the ground 
state to ground state transition. 

The present paper is organized as follows. In Section 2 we develop a model for the density- 
dependent nucleon-nucleon (NN) interaction in-medium. For the density-independent part, we 
use model space/renormalization group techniques to construct a low-momentum two-body 
nuclear interaction based on the N3L0 chiral NN potential. We then present analytic ex- 
pressions for the six unique density-dependent contributions to the in-medium NN interaction 
derived from the leading-order chiral 3NF at one-loop order. As suggested in ref. [14], the 
two low-energy constants, Co and Ce, of the medium- and short-range chiral 3NF are chosen 
to reproduce the binding energies of light nuclei for a given low-momentum scale Aio^-k- In 
Section 3 we present the results of our shell model calculation, and in Section 4 we give a brief 
summary and outlook. 

2 In-medium nuclear interaction 

2.1 Free-space low-momentum nucleon-nucleon interaction 

For many years there has been much effort devoted to the construction of realistic models 
of the nuclear interaction, which are important not only for understanding the properties of 
normal nuclei but also for constraining the structure of dense stellar objects such as neutron 
stars. The underlying assumption of these models is that the nuclear force arises primarily 
from the exchange of various mesons, though the choice of which mesons to include and how 
to account for the strong short-distance repulsion in the NN interaction are largely model- 
dependent. Nevertheless, various models containing between 30 to 40 free parameters (typically 
meson masses, coupling constants, and form factor cutoffs) are able to reproduce all of the 
experimental pp and pn scattering phase shifts below a laboratory energy of 350 MeV, as well 
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as the properties of the deuteron, with a per degree of freedom of about 1. 

Recently a program has been developed [15, 16] to understand the scale-dependence of the 
NN interaction from the point of view of effective field theory and the renormalization group. 
Since the short-distance behavior of the NN interaction is not constrained experimentally, 
renomalization group techniques have been used in refs. [15, 16] to evolve these (bare) NN 
interactions down to the scale at which our experimental information stops, that is, around a 
center of mass momentum p ~ Aio^.k ~ 2.1 fm~^. These low-momentum interactions, V^iow-fc; 
are phase-shift equivalent to the underlying bare interaction up to a predefined cutoff scale 
Aiow-k- As the decimation scale is reduced to Aiow-k — 2.1 fm~^ ~ 400 MeV, all of these low- 
momentum interactions flow to a nearly universal interaction. The method for constructing 
such an interaction is explained below. 

One begins with the half-on-shell T-matrix^ for free space scattering, 

T{p ,p) ^Vnn{p ,P) + -V ^ ^ q dq, (1 

71 Jo p^ — q^ 

from which one defines a low-momentum half-on-shell T-matrix by 

71ow-k(p ,P) = Mow-klP ,P) + -V / q dq, (2) 

TT Jo p — q 

where V denotes the principal value and the cutoff Aiow-k is a priori arbitrary. To preserve 
phase shifts, ta.n6{p) = —pT{p,p), one requires that these two T-matrices be equal for relative 
momenta less than the cutoff Aiow-k, 

Tiov,-k{p',p) = T{p',p) , for p',p < Aio„_k , (3) 

a condition which defines the low-momentum potential l^ow-fc- It can be shown [15] that a 
solution to these equations is given by the Kuo-Lee-Ratcliff folded diagram effective interaction 
[17, 18]. There are several schemes [19, 20] available for accurately computing Viow-k, and each 
scheme preserves the deuteron properties. Under this (scale) decimation procedure, all high- 
precision NN potentials flow, as Aiow-k — >^ 2.1 fm~^, to a nearly unique low-momentum potential 
Viow-k- In this study we consistently employ the N3L0 chiral NN potential in deriving Viow-k- 

Given this model of the nuclear interaction in free space, the question remains how to extend 
the description to a nuclear medium with densities close to that of saturated nuclear matter. It 
is well known that such two-body interactions alone are unable to reproduce simultaneously the 
known saturation energy and density of symmetric nuclear matter. The traditional approach 
is to include a three-nucleon force, but for many of the traditional NN interaction models it 
is difficult to systematically construct a 3NF that is consistent with the underlying two-body 
interaction. However, by exploiting the separation of scales in the framework of chiral effective 
field theory, a systematic and consistent construction of two-, three-, and four-nucleon forces has 
become possible (for a recent review, see ref. [21]). The key element there is a power-counting 
scheme which orders the contributions in powers of small external momenta over the chiral 
symmetry breaking scale. Long-range effects from multi-pion exchanges between nucleons are 
treated explicitly (and calculated within chiral perturbation theory), while the short-distance 
dynamics is encoded in nucleon contact terms. When applied to two- and few-nucleon problems, 
these chiral potentials are regulated by exponential functions [13, 21] with cutoffs ranging from 
500 to 700 MeV, in order to eliminate (unphysical) high-momentum components. 

The construction of decimated low-momentum three-body forces consistent with the two- 
body decimation for Viow-fe is currently a challenge. A common practice is to adjust the param- 
eters of the chiral 3NF so that the binding energies of ^H, ^He, and ^He are reproduced. Such 

^This real-valued quantity T{p',p) is often referred to as the JT-matrix. 
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a calculation was carried out in [14] using the Argonne Vis [22] low-momentum interaction, 
where in Table 1 we restate the values of the low-energy constants cd and ce determined in [14] 
as a function of Aiow-k- For a cutoff scale of Aiow-k = 2.3 fm~^, Viow-fc is only weakly model 
dependent. Therefore, we expect the values of cd and ce derived for the Argonne vis potential 
at Aiow_k = 2.3 fm^^ to be applicable for the low-momentum N3L0 chiral potential employed 
in the present paper. 



Alow-k 


Cd 


Ce 


2.1 fm-i 


-2.062 


-0.625 


2.3 fm"^ 


-2.785 


-0.822 



Table 1: The values of the low-energy constants cd and c^; of the chiral three- nuclcon interaction 
fit [14] to the binding energies of A = 3,4 nuclei for different values of the momentum cutoff 
Aiow-k- The interaction Viow-fc is derived from the Argonne Vig potential. 



2.2 In- medium nucleon-nucleon interaction 

In this section we derive from the leading-order chiral three-nucleon interaction an effective 
density-dependent in-medium NN interaction. As shown in ref. [23] the chiral three-nucleon 
interaction plays an essential role in obtaining the saturation of nuclear matter when using the 
universal low-momentum NN potential V^ow-k in Hartree-Fock calculations. The parameters 
of its two-pion exchange component, Ci = — O.SlGeV"^, C3 = — 3.2GeV'~^, C4 = 5.4GeV~^, 
are well determined from fits to low-energy NN phase shifts and mixing angles [13]. If one 
restricts the two large coefficients 03,4 to their (dominant) A(1232)-resonance contributions, 
the celebrated thrcc-nuclcon force of Fujita and Miyazawa [24] is actually recovered. The 
coupling parameters associated with the one-pion exchange component (cd) and the short- 
range contact-term (ce) of the chiral 3N-interaction have also been adjusted in refs. [14, 25] to 
binding energies of three- and four-nucleon systems (^H, ^He, and ^He). 

We are considering the (on-shcll) scattering of two nucleons in the nuclear medium in the 
ccntcr-of-mass frame, Ni{p) + N2{—p) Ni{p + q) + N2{—p — q), i.e. the total momentum of 
the two-nucleon system is zero in the nuclear matter rest frame before and after the scattering. 
The magnitude of the in- and out-going nucleon momenta isp = \p\ = \p + q\ , and ^=1^1 gives 
the magnitude of the momentum transfer. Having discussed the kinematics, we reproduce first 
(for the purpose of comparison) the expression for the (bare) l7r-exchange: 



. .(1^) _ 9a^n ^ ^<?i-go^2-g f.s 



Here, qa = 1-3 denotes the nucleon axial-vector coupling constant, f^^ = 92.4 MeV is the weak 
pion decay constant, and = 135 MeV is the (neutral) pion mass. Furthermore, ai,2 and ri,2 
are the usual spin and isospin operators of the two nucleons. Note that we have included an 
additional factor of Mjv/47r in Vnn in order to be consistent with the conventions chosen for 

V\ow-k- 

We start with those contributions to the in-medium NN-interaction V}^^ that arc gen- 
erated by the 27r-exchange component of the chiral three-nucleon force. The three different 
topologies for (non- vanishing) one-loop diagrams are shown in Fig. 1. The short double- line 
on a nucleon propagator symbolizes the filled Fermi sea of nucleons, i.e. the medium insertion 



4 




Figure 1: In-medium NN interaction generated by the two-pion exchange component (~ 01,3,4) 
of the chiral three-nucleon interaction. The short double-hne symbohzes the filled Fermi sea of 
nucleons, i.e. the medium insertion —27t6{Iq) 6{kf — |/ 1) in the in-medium nucleon propagator. 
Refiected diagrams are not shown. 



—2'k5{Iq) 9{kf — |/ 1) in the in-medium nucleon propagator. In effect, the medium insertion 
sums up hole propagation and the absence of particle propagation below the Fermi surface 
| / 1 < kf. The left diagram in Fig. 1 represents a Ivr-exchange with a Pauli blocked in-medium 
pion self-energy and the corresponding contribution to V^'^'^ reads: 



V, 



med,l _ 9a^NP 



NN 



TI-T2 



(Ti- qa2- q 

[ml + g2)2 



{2ciml + csq^) . 



(5) 



The Fermi momentum kf is related to the nucleon density in the usual way, p = 2fcj/37r^. 
Since ci,3 < 0, this term corresponds to an enhancement of the bare Ivr-exchange. It can be 
interpreted in terms of the reduced (spatial) pion decay constant, /*^^ = + 2c3p, in the 
nuclear medium that replaces in the denominator of eq.(5). The second diagram in Fig. 1 
includes vertex corrections to the Ivr-exchange caused by Pauli blocking in the nuclear medium. 
The corresponding contribution to the in-medium NN-interaction has the form: 



..med,2 _ 9aMn ^ ^ ^1 " ? ^ 2 " Q 
^NN — Qo 3/4''"! ■ ''"2 2~~i 2~ 

32vr'^/^ ml + q^ 

-(C3 + C4) [q\To{p) + 2T^{p) + T,{p)) + AT^ip) 



+ 4C4 



(6) 



i2l-l. 



Here, we have introduced the functions Tj{p) which result from Fermi sphere integrals over a 
(static) pion-propagator [ml + {I + p 

kf +p 



^o{p) = kf 
ri(p) 



m^ 



arctan ■ 



-|- arctan 



m^ 



m^ . 



ml + k'l 



4p 



J^^^ml + {kf+pf 



f 



-Up)- 



16p^ 



1.3 
9 



T2{p) = '^ + ^ikj-ml- p')Toip) + ^(m^ + A;J - p'')T,ip) 



ml + {kf+pY ml + {kf-p)" 
1 



ml + {kf — pY ' 

ml + {kf + pf 
ml + {kf - p)' 



In 



k'^ ml + k'f + p 



2p2 



-Up) - 



ml + kj + 



Up)- 



(7) 

, (8) 

(9) 

(10) 



By analyzing the momentum and density dependent factor in eq. (6) relative to V^jvtv'' ) finds 
that this contribution corresponds to a reduction of the Ivr-exchange in the nuclear medium. 
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Approximately, this feature can be interpreted in terms of a reduced nucleon axial- vector con- 
stant g*A{kf). 

The right diagram in Fig. 1 represents Pauh blocking effects on chiral 27r-exchange. Evaluat- 
ing it together with the reflected diagram one finds the following contribution to the in-medium 
NN-interaction: 



V, 



med,3 
NN 



N 



-C3 [9>k) - 12(2m2 + g2)ro(p) - Qq^T^{p) + 3(2m^ + q^fGo{p, q) 
-F4c4 fi • T2 {Bi • (72 - (Ti • g (72 • q)G2{p, q) 
-(3c3 + C4fi • 7^2) + <?2) • (g X p) [2ro(p) + 2ri(p) - {2ml + 
X {Go{p, q) + 2G,{p, q))] - Ucrml i{a^ + B^) ■ {q X p) [Go{p, q) + 2G,{p, q) 
+Aci fi • f2 ai • (g X p) a2-{qx p) \Gq{p, q) + 4Gi(p, q) + 4G3(p, g) 
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One observes that in comparison to the analogous 27r-exchange interaction in vacuum (see 
section 4.2 in rcf. [26]) the Pauli blocking in the nuclear medium has generated additional spin- 
orbit terms, i((7i + (J2) ■ {qx p), and quadratic spin-orbit terms, Bi ■ {q x p) ^2 ■ {q x p) , written 
in the last three lines of eq.(ll). The density dependent spin-orbit terms (scaling with 03^4) 
in the in-medium NN-interaction demonstrate clearly and explicitly the mechanism of 

three-body induced spin-orbit forces proposed long ago by Fujita and Miyazawa [24]. 

The functions Gj{p,q) appearing in eq.(ll) result from Fermi sphere integrals over the 
product of two different pion-propagators. Performing the angular integrations analytically 
one arrives at: 



0,*,** 



dl 



^^Mp) 



;i2) 



with the abbreviation A{p) = [m^ + [l+pY] [ml + [l—pY]- The other functions with j = 1, 2, 3 
are obtained by solving a system of linear equations: 



Giip,q) 

Gu{p,q) 
G2{p,q) 
Gz{p,q) 



ro{p)-{ml+p')Go{p,q)-G.{p,q) 
4p2 _ q2 ' 

3r2(p) +P^r3(p) - {ml+p')G4p, q) - G.4p, q) 
4p2 — g2 ' 

{ml + p^)Gi{p, q) + a(p, q) + Gu{p, q) , 

iri(p) - 2{ml + p2)G'i(p, q) - 2Gu{P: q) - G^p, q) 



(13) 

(14) 
(15) 
(16) 



_ g2 

In this chain of equations the functions indexed with an asterisk play only an auxiliary role for 
the construction of Gi,2,3(p, Q')- We note that all functions Gj{p,q) are non-singular ai q = 2p 
(corresponding to scattering in backward direction). 

Next, we come to the Ivr-exchange component of the chiral three- nucleon interaction pro- 
portional to the parameter where c^, ^ —2 for a scale of ~ 0.7 GeV [14]. The filled 
black square in the first and second diagram of Fig. 2 symbolizes the corresponding two-nucleon 
one-pion contact interaction. By closing a nucleon line at the contact vertex, one obtains a 
vertex correction (linear in the density p) to the Ivr-exchange: 

qaMnCdP ^ ^ Bi- qa2- q 



med,4 
NN 



327r/,4A, 



Tl • T2 



mi 



(17) 



6 




(6) 



Figure 2: In- medium NN interaction generated by the one-pion exchange 
range component (~ ce) of the chiral three-nucleon interaction. 



cd) and short- 



Since Cf) is negative, this term reduces again the bare Ivr-exchange, roughly by about 16% 
at normal nuclear matter density po = 0.16 fm~^. The second diagram in Fig. 2 includes 
Pauli blocked (pionic) vertex corrections to the short-range NN interaction. The corresponding 
contribution to the density dependent in-medium NN interaction reads: 



med,5 
NN 



ai ■ CT2 2p- 



cTi ■ g 0-2 ■ g ( 1 - 



— ^^i ■ (g X p)a2 ■ (g X p) 



ro(p) + 2ri(p) + r3(p) 



(18) 



where we have used an identity ioi Bi ■ p (72 ■ p + (Ti ■ {p + g) a2 ■ {p + q) = [■ ■ ]■ The ellipses 
stands for the combination of spin operators written in the square bracket of eq.(18). 

Finally, there is the short-range component of the chiral 3N interaction, represented by 
a three-nucleon contact- vertex proportional to ce/^x- -^y closing one nucleon line (see right 
diagram in Fig. 2) one obtains the following contribution to the in-medium NN-interaction: 



V, 



med,6 3MnCeP 



NN 



(19) 



2A;^/37r^ and is independent of spin, isospin and 



which simply grows linearly in density p 
nucleon momenta.^ 

For implementation of the in-medium NN-interaction into nuclear structure calculations, 
the matrix elements of in the LSJ basis are needed. These are obtained by setting 

Ti-T2 = 4T — 3, with T = 0, 1 the total isospin, and g = P\J2{1 — z) and performing a projection 

^\ dz Pl{z) . . . with Legendre polynomials (for details see section 3 in ref. [26]). The resulting 
diagonal spin-singlet {S = 0, L = J) and diagonal spin-triplet {S = 1, L = J— I, J, J+l) matrix 
elements as well as the off-diagonal mixing matrix elements {S = 1, L = J — 1, L' = J + 1) are 
then functions of the momentum p and kf (or equivalently the nucleon density p). Note that 
the off-diagonal mixing matrix elements arise exclusively from the tensor operator ai ■ qa2 ■ q 
and the quadratic spin-orbit operator ai ■ (g x p ) cr2 ■ (g x p). 

The restriction to on-shell amplitudes greatly simplifies the calculation of the medium- 
dependent corrections to the nuclear interaction. However, in calculating the shell model matrix 
elements one must know the components of the interactions also off-shell. That is, one needs 
to consider Ni{p) + N2{—p) Ni{p + g ) + N2{—p — g), where p = \p \ and p' = \p + q\ ^ P- 



■^In order to facilitate the computation of symmetry factors and spin and isospin traces, we have modeled 
(for that purpose) the three-nucleon contact-interaction by heavy isoscalar boson exchanges. 
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Figure 3: Two different off-sfiell extrapolations for the density dependent V^^'^ at p — po com- 
pared to the off- shell dependence of the low- momentum one-pion-exchangc in the ^5*0 partial- 
wave. Matrix elements are shifted and scaled for comparison purposes (see text for details). 



Rather than calculate the full off-shell in-medium NN interaction, which would give rise to 
new spin- dependent operators beyond those present in free-space two-nucleon scattering, we 
look at two prescriptions for extending the on-shell amplitudes off-shell. First, we consider 
the symmetric extrapolation obtained by replacing in the above on-shell expressions with 
|(p^ + p'^)- Second, we consider the asymmetric extrapolation obtained by simply introducing 
no explicit dependence on p'^. As we discuss later, the largest ofT-shcll dependence comes from 
V^^''^. In Fig. 3 we compare the two off-shell extrapolations for V^'^''^ at po = 0.16 fm~^ in the 
^Sq partial wave to the off-shell dependence of the one-pion-exchange (OPE) low-momentum 
interaction, where Aiow-k = 2.1 fm~^. We have fixed the value of p = 0.04 fm""^, and for 
ease of comparison among the three sets of off-diagonal matrix elements, we have changed the 
overall sign of the OPE interaction, shifted all three to the origin, and scaled the two medium- 
dependent interactions so that they coincide with the low-momentum OPE at p' = 1.6 fm~^. 
We find that there is excellent agreement between the symmetric off-shell extrapolation and 
the low-momentum OPE for momenta up to nearly p' — 1.75 fm~^. Therefore, throughout the 
remainder of this study we employ only the symmetric off-shell extrapolation. This completes 
our construction of the density-dependent in-medium NN interaction arising from the leading- 
order chiral three- nucleon force. 
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3 Shell model calculation 



3.1 Formalism 

Nowadays there are several computationally intense ab initio approaches to solve the nuclear 
many-body problem for finite nuclei. These include Green function Monte Carlo techniques 
[27], the no-core shell model [28], and coupled cluster methods [29, 30]. However, A = 14 nuclei 
are currently still too complicated for any of these methods to treat exactly. Therefore, in 
this study we employ the standard shell model, which restricts the allowable configurations to 
Ohu but includes the mixing of more complicated configurations through perturbation theory. 
Indeed, the shell model is expected to provide a good description of light nuclei close to a 
doubly-magic nucleus (e.g., ^^O in the present study). 

We describe the states of ^''C and ^^N as consisting of two Op holes in a closed ^^O core. 
The harmonic oscillator parameter is chosen to be hu = 14 MeV, which yields good agreement 
with the experimental charge distribution of ^^O. The second parameter in the model is the 
energy splitting between the two Op orbitals, e = e(pi/2) — 6(^3/2) = 6.4 MeV, which we have 
taken from the experimental excitation energy of the first 3/2^ state of ^^N located 6.4 MeV 
above the l/2~ ground state. 




Figure 4: Diagrams contributing to the effective interaction V^s in our present calculation. 
Wavy lines represent the density-dependent in-medium nuclear interaction caculated in Section 
2. 

To obtain the ground state energies and wavefunctions we construct the shell model effective 
interaction following the formalism explained in ref. [18]. The full nuclear many-body problem 

= K^n, (20) 

which presently cannot be solved for mass number A = 14, is replaced with a model space 
problem 

HeQ~X_m -l^mX,m- (21) 

In this equation 

H = Ho + V and = Hq + Ks, (22) 

where Hq is equal to the sum of the kinetic energy and single-particle harmonic oscillator 
potential. En = En{A = 14) — Eq{A = 16, core), and V denotes the input NN interaction 
(in this case, it is provided by the density-dependent potential Viow-k + Vnn^)- The effective 
interaction Ves is obtained from the folded- diagram formalism detailed in ref. [31]. In the 
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present study, the Q-box [31] is calculated using hole- hole irreducible diagrams of first- and 
second-order in Viow-k + shown in Fig. 4. To distinguish these many-body particle- 

hole diagrams from the pion-exchange diagrams contributing to the in-medium NN interaction 
calculated in Section 2, we symbolize the input interaction with a wavy line. In previous 
work [31], the strong short-distance repulsion in the NN S'-wave interaction was mitigated by 
constructing the G-matrix, thereby yielding an effective interaction suitable for perturbation 
theory techniques. Our present use of low-momentum interactions achieves the same purpose 
and obviates the construction of the G- matrix. Previous studies [23, 32, 33] have found that in 
general low-momentum interactions are suitable for perturbative calculations. In all of these 
references satisfactory converged results were obtained including terms only up to second order 
in Viow-fc- Moreover, low- momentum interactions are energy independent and do not suffer 
from possible double counting problems that can arise with G-matrix interactions. 
We denote the single-hole states by 



11) = 



Os 



-/,),|2)= 0p-/,},\3)^ Op-l,}, 



1 \ 



-1 \ 



(23) 



and the two-hole states coupled to good angular momentum J, isospin T, and parity tt are 
denoted by \aP; J^T), where a, /3 = 2 or 3. The matrices to diagonahze are 



(a/?;l+0|Veff|75;l+0) 




-e 







-2e 



(24) 



for the J'^ — 1+, T — states, and a similar 2x2 matrix for the J'^ 
ground states for -"^^C and ^^N are respectively written as^ 



ipi ^ a 22; 0+1) + 6 33; 0+1 



0+, T = 1 states. The 



X 



22; 1+0^ + y 23; 1+0^ + z 33; 1+0^ 
and the reduced Gamow-Teller matrix element Mqt is evaluated to be 

Mgt = E MH^Mkmi) = 4f (-2V5ax + 2V2ay + Aby + bz) . 
k v6 ^ 



(25) 



(26) 



The half-life T1/2 is inversely proportional to the square of the Gamow-Teller matrix element, 
and for the ^''C decay one has 



27r3^^1n2 



1/2 



f{Z,Eo) mlc^Gl gf\MGT 



2 ' 



(27) 



where Eq — 156 keV is the maximum kinetic energy of the emitted electron, f(Z,EQ) is the 
Fermi integral over phase space, Gy — Gp cosd^ is the weak vector coupling constant, and g'^ 
is the in-medium axial vector coupling constant. Although the small end-point energy Eq of 
the -"^^C beta decay suppresses the transition by a factor of several hundred, this is not nearly 
sufficient to account for the observed lifetime. In order to obtain the anomalously long half-life 
of Ti/2 ~ 5730 years, the GT matrix element must be anomalously small, Mgt — 2 x 10""^. 



''We note that the jj-coupling scheme employed in the present paper differs from the previous density- 
dependent study [8] in which the two-particle states were LS'-coupled. 
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3.2 Results 



As a preliminary step, we first calculate the ^^C and ^^N wavef unctions, as well as the Gamow- 
Tcllcr transition matrix element, using the bare Viow-fe derived from the Idaho N3L0 potential 
at a cutoff scale Aiow-k = 2.1 fm~^. We find 



J'^ = 0+, T = 1 
jTT ^ 0+, T = 1 



0.395 
0.919 



22; 0+,l) + 0.919 33;0+,l) , 
22; 0+,l) - 0.395 33;0+,l) , 



(28) 



for the two 0+ states of ^^C with energy splitting A£^(0"'") — 12.9 MeV. For the two lowest 1+ 
states of ^''N we find 



i+,r = 

1+,T = 



0.137 
0.360 



22; 0+, l) - 0.676 33; 0+, l) + 0.725 33; 0+, l) , 

22; 0+, l) - 0.670 33; 0+, l) - 0.649 33; 0+, l) , (29) 



where the first excited state lies AE{1~^) = 2.77 MeV above the ground state. The second 
excited J'^ = 1~^,T = state lies nearly 16 MeV higher in energy than the ground state and 
will therefore be neglected throughout. With these wavef unctions, the ground state to ground 
state Gamow-Teller transition matrix element is found to be 



Mgt = -0.88, 

which is much too large to describe the known hfetime of ^^C. 



(30) 



j^^o+,r^i (p = po/io) 





(22|VV|22) 


(2211/^^7^-^133) 


(33|V77^^|33) 


1 


0.277 


0.053 


0.240 


2 


-0.292 


-0.071 


-0.242 


3 


0.043 


-0.048 


0.077 


4 


-0.054 


-0.011 


-0.047 


5 


0.041 


0.027 


0.022 


6 


0.171 


0.121 


0.085 



r = 0+,T = l (p = Po) 





(2211/^^^,-^122) 


(22|\/^f^,<^|33) 


(33|\/^f^.<^|33) 


1 


2.76(3 


0.525 


2.395 


2 


-3.367 


-0.714 


-2.863 


3 


1.231 


0.004 


1.228 


4 


-0.545 


-0.107 


-0.469 


5 


0.458 


0.316 


0.235 


6 


1.707 


1.207 


0.853 



Table 2: Matrix elements (in units of MeV) between Op ^ states coupled to {J^, T) = (0+, 1) for 
the six density-dependent contributions to the in-medium NN interaction at the two densities 
Po/10 and po- 



We now discuss the role of the six density-dependent components V^^-^ derived in Section 
2 that arise from the lowest-order chiral 3NF. Indeed, since A — lA nuclei lie just below a 
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double shell closure, we expect the average density experienced by a valence p-shell nucleon to 
be close to that of saturated nuclear matter. To begin we consider only the matrix elements 
of these interactions within the Op~^ model space. These are shown in Tables 2 and 3 for each 
of the six V/vA^^'* at the two densities p — po/10 and p = pQ. The matrix elements connecting 
states outside the model space contribute at second order in the diagrammatic expansion of the 
effective shell- model interaction, V^g, and will be included later. In general, the largest matrix 



J- = 1+,T = (p = Po/10) 





(22|l^;f^'i|22) 


{2?>\V^^^\22) 


(23|\/^r|23) 


(23|\/^r|33) 


(22|\/^^^'i|33) 


(33|\/^^^|33) 


1 


0.299 


-0.051 


0.128 


0.143 


-0.121 


0.240 


2 


-0.257 


0.091 


-0.147 


-0.101 


0.151 


-0.242 


3 


0.056 


0.028 


-0.019 


-0.015 


-0.100 


0.182 


4 


-0.054 


0.012 


-0.027 


-0.023 


0.025 


-0.047 


5 


0.034 


-0.026 


0.012 


0.021 


-0.040 


0.034 


6 


0.103 


-0.107 


0.171 


-0.002 


-0.055 


0.089 





(22|y^^^|22) 


(23|y^j^d|22) 


(23|y^^^|23) 


(23|y^^^|33) 


(22|y^^^|33) 


(33|y^j^^|33) 


1 


2.994 


-0.505 


1.282 


1.427 


-1.209 


2.395 


2 


-3.193 


0.865 


-1.718 


-1.291 


1.589 


-2.863 


3 


0.634 


-0.464 


0.524 


-0.035 


-0.927 


1.255 


4 


-0.542 


0.123 


-0.273 


-0.230 


0.247 


-0.469 


5 


0.322 


-0.286 


0.319 


0.100 


-0.276 


0.303 


6 


1.027 


-1.074 


1.714 


-0.016 


-0.550 


0.885 



Table 3: Matrix elements (in units of MeV) between Op^^ states coupled to ( J'^, T) = (1+, 0) for 
the six density-dependent contributions to the NN interaction in-medium at the two densities 
Po/10 and po- 



elements are those arising from V^^'^ and V^vw^'^' which are respectively the in-medium pion 
self-energy and vertex correction resulting from the long-range chiral 3NF. However, we find 
that to a good approximation they cancel. The same is true for V^^'^ and V^^'^ , which are 
relatively small to begin with. Therefore, we would expect V^^*^'^ and V^'^'*^ to be the most 
important density-dependent corrections to the in-medium NN interaction for the observables 
considered here. 

These observations are made more concrete in the following. As a measure of the strength 
of the various components contributing to the density-dependent part of the NN interaction, 
we define the following quantity 

{al3J-T\V,o^.^\^6J-T) ' ^^^^ 

where p is the nuclear density. The sum over a, (3, 7, and 6 runs over all possible configurations 
with the allowed spin, isospin, and parity, and n is the number of such configurations (e.g., n = 3 
for the J'' = 0+,T = 1 states). We derive Viow-fc from the chiral N3L0 NN interaction and 
choose Aiow-k = 2.1 fm~^. A value of R- close to 1 indicates that the relevant matrix elements 
for a particular medium correction are on average equal in magnitude to the matrix elements 



n 
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2r^ < 1 1 1 < < < 1 1 < < < 1 1 < < < 1 1 < < < 1—1 1| — < — < — < — < — I — < — < — < — < — I — < — < — < — < — I — < — < — < — < — I — < — < — < — r 




Figure 5: The average strength of the three (a/3; 0+l| V^^'^ \'yS] 0^1) matrix elements for the 
six density-dependent contributions to the in-medium NN interaction relative to Vjow-k- The 
low-momentum interaction is constructed from the Idaho N3L0 potential for Aiow-k = 2.1 
fm~^. 



of Viow-k in the given spin-isospin state. We show on the left side of Figs. 5 and 6 the value of 
R for each of the six density- dependent components of the in-medium nuclear interaction for 
densities between p = and p = po- The plots on the right hand side of the two figures include 
also the ratio R for V'^^'^'^ + V^""^'^ and V'^^'^'^ + \/^«d,5_ ^^lat V^^"^'^ + V^^'^''^ as well as 

V^j^*^'^ + ^N^'^ largely cancel, so that the most important contributions from the chiral 3NF 
are the Pauli-blocked 27i exchange and the three-nucleon contact interaction proportional to 

Ce- 

Before discussing the results of the full calculation shown at the end of this section, we first 
obtain a physical understanding of the role played by the two dominant components, y™'^^'^ and 
ymedfi^ from lowcst-order perturbation theory. Here the density-dependent components induce 
small changes in the ground state wavefunctions by mixing in the first excited states. This is 
justified on the basis that the relevant matrix elements of the density-dependent interactions 
are small compared to the lowest-order matrix elements from V^ow-k as seen in Figs. 5 and 6. 
We include explicitly the higher-order many-body diagrams shown in Fig. 4 at the end of this 
section. 

Given the ground state wavefunctions obtained from \^ow-k alone, as shown in eqs. (28) 
and (29), we see that in eq.(26) the first three terms contributing to the GT matrix element 
enter with the same sign. It is only the last contribution {bz) that reduces the strength of the 
transition matrix element. In order to suppress the GT transition, the density-dependent 
corrections must therefore shift strength to the |33) components of the wavefunctions. A 
straightforward calculation using the matrix elements shown in Tables 2 and 3 for p = po 
gives for the un-normalized perturbed wavefunctions arising from V^^'^ 

4^)(0+,l) = 4°)(0+,l)- 0.0003^^(0^,1), 

4^)(l+,0) = 4°)(l+,0) + 0.002V^f)(l+,0), (32) 



and similarly from V^^'^ 



rf)(0+,l) = 4°)(0+,l)- 0.088^^^(0+, 1) 
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0.4 0.6 
p/Po 



0.4 0.6 

p/Po 



Figure 6: The average strength of the six (a/?; l+Ol V^^'^ \'yS; 1~^0) matrix elements for the 
six density-dependent contributions to the in-medium NN interaction relative to Vjow-k- The 
low-momentum interaction is constructed from the Idaho N3L0 potential for Aiow-k = 2.1 
fm~^. 



4'\l\0) = V^r(1^0)-0.313<^(l+,0). 



,(0) 



(0), 



(33) 



This leading perturbative calculation shows the following features. First, the effects due to 

and second, one expects a reduction in the GT transition 
is to increase the strength of the |33) components of both 



V^^'^ dominate those from V^^y > 
strength, since the effect of V^^'^ 
ground state wavefunctions. 

We now consider the results of the full calculation, including all density- dependent contri- 
butions to second-order in the shell model effective interaction shown in Fig. 4. In Table 4 we 
show the results of our calculations for the expansion coefficients of the ground state wavefunc- 
tions of ^^C and ^^N, as well as the reduced GT matrix element, as a function of the nuclear 
density up to p = 1.25po- We see that the main effect of including the density-dependent 



p/po 


a 


b 


X 


y 


z 


Mgt 


0.00 


0.40 


0.92 


0.14 


-0.68 


0.72 


-0.877 


0.25 


0.37 


0.93 


0.13 


-0.67 


0.73 


-0.833 


0.50 


0.34 


0.94 


0.11 


-0.63 


0.77 


-0.684 


0.75 


0.30 


0.95 


0.09 


-0.57 


0.82 


-0.488 


1.00 


0.25 


0.97 


0.07 


-0.49 


0.87 


-0.267 


1.25 


0.19 


0.98 


0.05 


-0.41 


0.91 


-0.045 



Table 4: The coefficients of the jj-coupled wavefunctions defined in eq. (25) and the associated 
reduced GT matrix element as a function of the nuclear density p. 



modifications of the NN interaction is to strongly suppress the GT transition, as we expected 
from our previous perturbative analysis. However, in this calculation with theoretical errors 
neglected, it appears that only at densities above that of saturated nuclear matter would the 
suppression be strong enough to reproduce the experimentally observed half-life. 
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The Gamow- Teller transition strength, B{GT), is related to the reduced GT transition 
matrix element by 

B{GT) = {gX)'^j^\MGT\' ^ iMcTr, (34) 

where Jj = and we have used the approximation that the in-medium axial vector coupling 
constant (7^ ~ 1. In general, the effective Gamow- Teller operator has additional terms in a 
nuclear medium [34, 35]: 

^GT,cff = QlaL + g*A^ + gpA [>2, (35) 

where L is the orbital angular momentum operator and Y2 denotes the rank-2 spherical harmonic 
y2,m- However, it is well known from theoretical calculations [36] of the effective GT operator 
in-medium, as well as beta decay calculations [37] performed with phenomenological shell model 
effective inter ations, that g^A and qpa are almost negligible and that in light nuclei g\ is smaller 
by 15-20% from its free space value of qa = 1-27 (measured in neutron beta decay). In general, 
one must calculate both the effective interaction V^s and the effective Gamow- Teller operator 
CGT,eflf- In our calculations we assume a 20% reduction of qa in medium and therefore set 
9*A = 1-0. 

In Fig. 7 we separate the effects of the six different components of V^jv* on the GT strength. 
Including only the sum of V^^'^, V^at^'^, V)viv^'^i and V^^'^'^ together with the low-momentum 
interaction Viow-k) there is very little density dependence in the calculated GT strength. These 
four density- dependent components are exactly the terms that are expected to approximately 
cancel at leading order in the effective shell model interaction. Including the additional V^^'^, 
we find that there is a mild increase in the GT strength. This is in accordance with our 
predictions based on the first-order perturbative calculation of the ^^C and ^^N ground state 
wavefunctions shown in cq.(32). However, given the very small mixing of the excited state 
wavefunctions at first-order, the second-order diagrams contributing to the shell model effective 
interaction are equally important. Finally, by including V^^^'^ in addition to all other terms, 
we find that the GT strength is strongly suppressed. From eq.(33) our interpretation is that 
the 3N contact interaction strongly shifts strength to the |33) components of both the ^'^C and 
i^jvj ground state wavefunctions, thereby leading to the observed suppression. 

In Section 2 we suggested that the symmetric off-shell extrapolation of the density-dependent 
interaction is preferable to the asymmetric extrapolation, given the former's close agreement 
with the off-shell behavior of the OPE Viow-k- In Fig. 8 we study the effect of this choice on 
the Gamow- Teller strength. Regardless of the off-shell extrapolation method, we find a clear 
suppression of the B{GT) value for densities close to that of nuclear matter. Indeed, the most 
important component of the density- dependent NN interaction, V^vjv^'^) a contact interaction 
and is therefore independent of the off-shell extrapolation used. We conclude that the error as- 
sociated with the off-shell prescription is rather small, given that even the extreme asymmetric 
extrapolation is rather close to the results from the symmetric extrapolation. 

Given the strong dependence of the ^^C beta decay transition rate on the nuclear density, 
it is important to study how other Gamow- Teller strengths are effected by the inclusion of 
V^^. Recently, GT strengths from the ^''N ground state to excited states of ^^C have been 
determined from the experimental charge exchange reaction ^^N(d, ^He)^^C [11]. In Fig. 9 
we plot our results for the calculated B{GT) values as a function of the nuclear density p. 
From Fig. 9 one sees that the medium effects improve the agreement between our theoretical 
calculations and the experimental values for transitions from 0"*" and 1"*" states. The largest 
effect is clearly a suppression of the ground state to ground state transition for densities at 
and above that of nuclear matter. The other transition strengths are much less sensitive to the 
density dependence of the nuclear interaction. 
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Figure 7: The effects of the various density-dependent contributions to the in-medium NN 
interaction on the ground state to ground state B{GT) value. The legend denotes which of the 
V^^'^ are included together with Viow-fc- 



Although we have chosen a low- momentum cutoff scale of Aiow-k = 2.1 fm~^ in the above 
calculations, the value of c^; (and consequently the strength of the important V}!^^'^ term) de- 
pends sensitively on Aiow-k as seen in Table 1. We therefore calculate the beta decay transition 
strength using also a cutoff scale of Aiow-k = 2.3 fm~^, which incorporates some components 
of the bare NN interaction unconstrained by elastic NN scattering data but which provides an 
estimate for the error associated with the chosen momentum cutoff. The results for densities 
up to 1.25po are shown in Fig. 10 together with the calculations using Aiow~k = 2.1 fm^^. We 
suggest that this constitutes the most important source of error in our calculation. Additional 
errors include the following: deviations of the effective axial coupling constant g\ from 1.0, the 
choice of the off-shell extrapolation for the in-medium NN scattering amplitudes, and finally, 
errors associated with our use of the shell model at second-order in perturbation theory. In 
earlier sections of the paper we have estimated errors associated with the first two approxima- 
tions, but the last of the three is difficult to quantify. Nevertheless, we expect such errors to be 
minor, given that previous shell model calculations (see [38] and references therein) using Vfow-fe 
up to second order in perturbation theory have been able to describe very well the properties 
of nuclei with a small number of valence nucleons above a closed shell. 

Since the theoretical value of B{GT) is particularly sensitive to the nuclear density, let us 
estimate the range over which we expect our calculations to be valid. In Fig. 11 we plot twice 
the charge distribution of ^"^N obtained from electron scattering experiments [39, 40] fit to the 
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Figure 8: Dependence of the ground state to ground state B{GT) value on the off-shell ex- 
trapolation of the density- dependent V^ff interaction. Both the symmetric and asymmetric 
off-shell extrapolation are shown. 



harmonic oscillator density distribution 

n(r)oc l^ + b"^ (36) 

where for ^^N, h = 1.291 and d = 1.740 fm. We compare this density distribution to the 
square of the Op-shell radial wavefunction used in our calculation. This wavefunction peaks at 
a nuclear density of approximately p = 0.85po- We have therefore shaded the region of Fig. 10 
corresponding to this density ± approximately 10%. 



4 Summary and outlook 

We have studied the effects of the leading-order chiral three-nucleon force on the beta decay 
lifetime of ^^C, and more generally on the Gamow- Teller transition strengths from the ground 
state of ^^N to the low-lying states of ^^C. Our results indicate that the density-dependent 
in-medium NN interaction V§^^ derived from the chiral 3N force has a strong effect on the 
ground state to ground state transition. In contrast, the GT strengths from the ground state 
of ^^N to the excited states of ^^C exhibit only a small density dependence. These results are 
consistent with the calculations presented in [8]. We find that the GT strength is particularly 
sensitive to the genuine short-range component of the chiral 3NF (which is almost completely 
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Figure 9: The B{GT) values for transitions from the low-lying states of ^^C to the ground state 
of ^^N as a function of the nuclear density. The experimental values are from [11]. Note that 
there are three experimental low lying 2"^ states compared to two theoretical 2"^ states in the 
Op~^ configuration. 



responsible for driving the suppression) as well as the low-momentum decimation scale Aiow-k- 
In general, by fitting the binding energies of A = 3,4 nuclei, one can derive a constraint curve 
relating the parameters cd and ce [41]. Additional information is then required to fix the point 
on this curve. Given the sensitivity of the ^^C lifetime on the parameter c^, we suggest that 
this decay can serve as a useful constraint on the two low-energy constants and ce once ah 
initio many-body calculations are able to study this problem more accurately. 
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0.25 0.5 0.75 1 1.25 



Figure 10: The uncertainty in the calculated value of B{GT) obtained by varying the low 
momentum cutoff Aiow-k between 2.1 and 2.3 fm~^. The shaded region corresponds to nuclear 
densities close to that experienced by valence Op-shell nucleons in ^^C. 
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Figure 11: Twice the charge distribution of ^^N taken from [39, 40] together with the square of 
the radial Op-shell wavefunctions. 
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